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Abstract 

Based on the Relativistic continuum Hartree-Bogoliubov (RCHB) theory, the 
pseudo-spin approximation in exotic nuclei is investigated in Zr and Sn isotopes 
from the proton drip line to the neutron drip line. The quality of the pseudo-spin 
approximation is shown to be connected with the competition between the centrifu¬ 
gal barrier (CB) and the pseudo-spin orbital potential ( PSOP ). The PSOP depends 
on the derivative of the difference between the scalar and vector potentials dV/dr. If 
dV/dr = 0, the pseudo-spin symmetry is exact. The pseudo-spin symmetry is found 
to be a good approximation for normal nuclei and to become much better for exotic 
nuclei with highly diffuse potential, which have dV/dr ~ 0. The energy splitting of 
the pseudo-spin partners is smaller for orbitals near the Fermi surface ( even in the 
continuum ) than the deeply bound orbitals. The lower components of the Dirac 
wave functions for the pseudo-spin partners are very similar and almost equal in 
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magnitude. 


PACS numbers : 21.10.Hw, 21.60.-n, 21.10.PC, 21.60.Jz, 27.60.j 


1. Introduction 


In a recent letter!, by relating the pseudo-spin symmetry back to the Dirac equation through 
the framework of Relativistic Continuum Hartree-Bogoliubov (RCHB) theory!, the pseudo¬ 
spin approximation in real nuclei was discussed. From the Dirac equation, the mechanism 
behind the pseudo-spin symmetry was studied and the pseudo-spin symmetry was shown 
to be related with the competition between the centrifugal barrier (CB) and the pseudo¬ 
spin orbital potential ( PSOP ), which is mainly decided by the derivative of the difference 
between the scalar and vector potentials. With the scalar and vector potentials derived from 
a self-consistent RCHB calculation, the pseudo-spin symmetry and its energy dependence 
have been discussed!. Here we will extend our previous investigation! to exotic nuclei. 
The pseudo-spin symmetry approximation for exotic nuclei is investigated for Zr and Sn 
isotopes ranging from the proton drip line to the neutron drip line. The isospin and energy 
dependence of the pseudo-spin approximation are investigated in detail. 

The concept of pseudo-spin is based on the experimental observation that the single 
particle orbitals with j = 1 + 1/2 and j = (l + 2) — 1/2 lie very close in energy and can 
therefore be labeled as pseudo-spin doublets with quantum number h = n— 1, 1 = 1 — 1, 
and s — s — 1/2. This concept was originally found in spherical nuclei 30 years ago!’!, but 
later proved to be a good approximation in deformed nuclei as well!. It is shown that the 
pseudo-spin symmetry remains an important physical concept even in the case of triaxialityi 


Since the suggestion of the pseudo-spin symmetry, much efforts has been made to un¬ 
derstand its origin. Apart from the rather formal relabeling of quantum numbers, various 
proposals for an explicit transformation from the normal scheme to the pseudo-spin scheme 
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have been made in the last twenty years and several nuclear properties have been investigated 
in this scheme! 0. Based on the single particle Hamiltonian of the oscillator shell model the 
origin of pseudo-spin was proved to be connected with the special ratio in the strength of 
the spin-orbit and orbit-orbit interactionJiE and the unitary operator performing a trans¬ 
formation from normal spin to pseudo-spin space was discussed0~E3. However, it was not 
explained why this special ratio is allowed in nuclei. The relation between the pseudo-spin 
symmetry and the relativistic mean field ( RMF ) theoryEl was first noted in Ref.!, in 
which Bahri et al found that the RMF explains approximately the strengths of spin-orbit 
and orbit-orbit interactions in the non-relativistic calculations. In a recent paper Ginocchio 
took a step further and revealed that pseudo-orbital angular momentum is nothing but the 
“orbital angular momentum” of the lower component of the Dirac wave functiorE. He also 
built the connection between the pseudo-spin symmetry and the equality in the scalar and 
vector potentialEE. 

To understand to what extent it is broken in real nuclei, some investigation along this 
line has been done for square well potentials^ and for spherical solutions of the RMF 
equationE. By relating the pseudo-spin symmetry back to the Dirac equation through 
the framework of RCHB theory, the pseudo-spin approximation in real nuclei was shown to 
be connected with the competition between the centrifugal barrier (CB) and the pseudo¬ 
spin orbital potential ( PSOP ), which is mainly decided by the derivative of the difference 
between the scalar and vector potentials. With the scalar and vector potentials derived from 
a self-consistent Relativistic Hartree-Bogoliubov calculation, the pseudo-spin symmetry and 
its energy dependence have been discussed in Ref.0. 

The highly unstable nuclei with extreme proton and neutron ratio are now accessible 
with the help of the radioactive nuclear beam facilities. The physics connected with the 
extreme neutron richness in these nuclei and the low density in the tails of their distributions 
have attracted more and more attention not only in nuclear physics but also in other fields 
such as astrophysicsSE. New exciting discoveries have been made by exploring hitherto 
inaccessible regions in the nuclear chart, ft is very interesting to investigate the pseudo-spin 
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symmetry approximation both in normal and exotic nuclei. For this purpose, we will use 
the RCHB theory, which is the extension of the RMF and the Bogoliubov transformation in 
the coordinate representation, and provides not only a unified description of the mean field 
and pairing correlation but also the proper description for the continuum and its coupling 
with the bound statefE. As this theory takes into account the proper isospin dependence 
of the spin-orbit term, it is able to provide a good description of global experimental data 
not only for stable nuclei but also for exotic nuclei throughout the nuclear chart!!. It is very 
interesting to examine the pseudo-spin symmetry approximation in exotic nuclei, in which 
the mean field potentials are expected to be highly diffuse. 

Recently, by relating the pseudo-spin symmetry back to the Dirac equation through the 
framework of RCHB theory, the pseudo-spin approximation in real nuclei was discussed. The 
mechanism behind the pseudo-spin symmetry was studied and the pseudo-spin symmetry 
was shown to be connected with the competition between the CB and the PSOP, which is 
mainly decided by the derivative of the difference between the scalar and vector potentials. 
With the scalar and vector potentials derived from a self-consistent RHB calculation, the 
pseudo-spin symmetry and its energy dependence have been discussed^. Here we will extend 
the previous investigation to the case of exotic nuclei. The pseudo-spin splitting in Zr and 
Sn isotopes has been studied from the proton drip line to the neutron drip line. The energy 
splitting of the pseudo-spin partners, their energy and isospin dependence will be addressed. 
An outline of the RCHB formalism is briefly reviewed in Sec. II. In Sec. Ill, the Dirac 
equation and the formalism leading to the pseudo-spin symmetry is presented. The energy 
splitting of the pseudo-spin partners and its energy dependence are given in Sec. IV. The 
pseudo-spin orbital potential, which breaks the pseudo-spin symmetry will be studied in Sec. 
V. In Sec. VI, the wave-function of pseudo-spin partners will be studied. A brief summary 
is given in the last section. 
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2. An outline of RCHB Theory 


The RCHB theory is obtained by combining the RMF and the Bogoliubov transformation 



in the coordinate representation^], and its detailed formalism and numerical solution can be 
found in Ref.i and the references therein. The RCHB theory can give a fully self-consistent 
description of the chain of Lithium ranging from 6 Li to n Li. The halo in n Li 

has been successfully reproduced in this self-consistent picture and excellent agreement with 



recent experimental data is obtained. The contribution from the continuum has been taken 
into account and proved to be crucial to understand the halo in exotic nuclei. Based on 


the RCHB, a new phenomenon ’’Giant Halo” has been predicted. The ’’Giant Halo” is 


composed not only of one or two neutrons, as is the case in the halos in light p-shell nuclei, 



the shell structure are systematically studied with RCHB in Ref.t-1 where both the pairing 
and blocking effect have been treated self-consistently. Therefore the RCHB theory is very 
suitable for the examination of the pseudo-spin approximation in exotic nuclei. 

The basic ansatz of the RMF theory starts from a Lagrangian density by which nucle¬ 
ons are described as Dirac particles interacting via the exchange of various mesons and the 
photons. The mesons considered are the scalar sigma (cr), vector omega (u>) and iso-vector 
vector rho (p). The iso-vector vector rho (p) meson provides the necessary isospin asymme¬ 
try. The scalar sigma meson moves in the self-interacting field of cubic and quadratic terms 
with strengths g 2 and gs, respectively. The Lagrangian then consists of the free baryon and 
meson parts and the interaction part with minimal coupling, together with the nucleon mass 
M, and m a , g a , m u , g U) m p , g p the masses and coupling constants of the respective mesons: 


C = - M)-0 + \d p ad^a — U(a) — 

+ \™g>Pv$ L ~ \ F ^ F ^ 
-g- g^ipip - gpi’pTijj - e^/H/L 


( 1 ) 


The field tensors for the vector mesons are given as: 


5 


iv w = av-av, 

Rf w = d^ - d u f?‘ - g P{pP x ?), 

piii' _ QHJp _ Qv 2 \fi' 


( 2 ) 


For a realistic description of nuclear properties, a nonlinear self-coupling of the scalar mesons 
turns out to be crucialil: 


U(a) = -m^cr* + — a° + —a‘ 

\ j 2 3 4 

The classical variation principle gives the following equations of motion : 
[a-p + V v {f) + /3(M + V s (r))}^i = 
for the nucleon spinors and 


(3) 


(4) 


(—Act + U'(a)) 

9aPs 

(—A + mfyu)* 1 

= 9uf(f) 

(-A + mfj Y 

= 9 P ff) 

-A Aq(t) 

= ej£(r) 


(5) 


with U'(cr) = d a U(c r) and A = —d^d^ for the mesons, where 


Vv(r) = + g p pr +-e(l - t 3 )4, 


( 6 ) 


V s (r) = g<jCr(r) 

are the vector and scalar potentials respectively and the source terms for the mesons are 


Pa = Ei=1 AA 
Af) = EtiAYA 
A if) = EiLi AY^A 
jiff = Dii lAY^-rp-A, 


(7) 


where the summations are over the valence nucleons only, ft should be noted that as 
usual, the present approach neglects the contribution of negative energy states, i.e., no¬ 
sea approximation, which means that the vacuum is not polarized. The coupled equations 
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Eq.(f|) and Eq.(||) are nonlinear quantum field equations, and their exact solutions are very 
complicated. Thus the mean field approximation is generally used: i.e., the meson field 
operators in Eq.(|4|) are replaced by their expectation values, so that the nucleons move 
independently in the classical meson fields. The coupled equations are self-consistently 
solved by iteration. 

For spherical nuclei, i.e., the systems with rotational symmetry, the potential of the 
nucleon and the sources of meson fields depend only on the radial coordinate r. The spinor 
is characterized by the quantum numbers Z, j,m, and the isospin t = ±— for neutron and 
proton, respectively. The other quantum number is denoted by i. The Dirac spinor has the 
form: 


ift(r) = 




,G ‘— N 
(a-~r)Y‘J9,<p) 


Jf|r) 

y r 




( 8 ) 


where Yj m (9, (ft) are the spinor spherical harmonics and G l f{r) and F 9 {r) are the radial wave 
function for upper and lower components. They are normalized according to the relation: 


/Dr(|Gf(r)| 2 + |if(r)| 2 ) = l. 
Jo 


(9) 


The radial equation of spinor Eq. 


can be reduced as 


<*G?(r) = + ^)if (r) + (M + V s (r) + V v (r))GV (r) 

or r 

e,F‘>(r) = (+1- + ^)Gt(r) - (M + Vs(r) - Vv{r))FV(r), 


( 10 ) 


dr r 


where 


k = 


~(j + 1/2) for j = l + 1/2 


+ {j + 1/2) for j = 1- 1/2. 

The meson field equations become simply radial Laplace equations of the form: 


d 2 


dr 2 


2 9 2 \ , 

r^ + ? M 0 = ^ (r) ’ 


( 11 ) 


are the meson masses for (ft — a, cu,p and for photon ( m$ = 0 ). The source terms are: 
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s<l>( r ) = * 


( 12 ) 


-gaps - 92 <? 2 {r) 

Qujpv 

gpPsir) 

ep c {r) 


( 7 3 (T 3 (r) for the a field 
for the u held 
for the p held 
for the Coulomb held, 


4vr r 2 p s (r) = Ef =1 (|^(r)| 2 - |F,(r)| 2 ) 
4 tt r 2 p v (r) = Eti(\ G i( r )\ 2 + \ F i( r )\~) 


r 2 p 3 (r) = E^HG p (r)j 2 + |F p (r)| 2 ) - £" =1 (|G n (r)| 2 + |F„(r)| 2 ) 


(13) 


4nr 2 p c (r) = E^HG p (r)j 2 + |F p (r)| 2 ). 

The Laplace equation can be solved by using the Green function: 

roo 

0(r) = / r' 2 dr'G ^(r , r') s ^r’), 

Jo 

where for massive helds 

GJr,r') = ———( e - m ^ r - r '\ _ e - m ^ r+r 'l) 

2 rr' 

and for Coulomb held 


(14) 


(15) 


G>(r,r') = 


1/r for r > r' 


(16) 


1/r' for r < r'. 

The Eqs.(jlOj) and ([11]) conld be solved self-consistently in the usual RMF approximation. 
However, Eq.(H)j) does not contain the pairing interaction, as the classical meson helds are 
used in RMF. In order to have the pairing interaction, one has to quantize the meson 
helds which leads to a Hamiltonian with two-body interaction. Following the standard 
procedure of Bogoliubov transformation, a Dirac Hartree-Bogoliubov equation could be 
derived and then a unified description of the mean held and pairing correlation in nuclei 
could be achieved. For the details, see Ref.! and the references therein. The RHB equations 
are as following: 

^ / tllTT \ / lllTT \ 

(17) 


d 3 r' 


h- A A 
A -h + A 


0 U 

,0V, 


= E 


0V 

,0V, 
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where 


h(r, r') = [a- p + V v {r) + /3(M + V s (r))} 8(r, r ) 


(18) 


is the Dirac Hamiltonian and the Fock term has been neglected as is usually done in RMF. 
The pairing potential is : 


A kk'(ry) = - j d 3 ri j 

w 


(19) 


It is obtained from the one-meson exchange interaction V kk , kk i(ZP\ fyE]) in the pp-channel 
and the pairing tensor k = V*U T \ 


Kkk'(ry) =< \a k a k ’\ >= if) T 


The nuclear density is as following: 


p(r,r) = J29iij^v(r)*^v(^)- 

ilj 


( 20 ) 


( 21 ) 


As in Ref!, V used for the pairing potential in Eq. (|TD|) is either the density-dependent two- 
body force of zero range with the interaction strength Vo and the nuclear matter density 
Po- 


V(ri, r 2 ) = Vq8(ti - r 2 )y [1 - (Tier 2 ] (l - 

4 V Po 


( 22 ) 


or 


Gogny-type finite range force with the parameter /i, : , l-'F,, BHi and Mi (i = 1, 2) 


V(r 1 ,r 2 ) = e 

i= 1,2 


((T’ 1 -'r 2 )//ii ) 2 


(Wi + BiP a - H t P T - MiP a P 1 


(23) 


A Lagrange multiplier A is introduced to fix the particle number for the neutron and proton 
as N = Trp n and Z = Tr p p . 

In order to describe both continuum and bound states self-consistently, we use the RHB 
theory in coordinate representation, i.e., the Relativistic Continuum Hartree-Bogolyubov ( 
RCHB ) theoryi. It is then applicable to both exotic nuclei and normal nuclei. In Eq. (0. 
the eigenstates occur in pairs of opposite energies. When spherical symmetry is imposed on 
the solution of the RCHB equations, the wave function can be written as: 
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(24) 


tij = 


( .Gl j (r) \ 


( < j (r) \ 


F$ j ( 


\ 


-((j-f) 


y l jm^A)Xt{t)^V = 


Fy j \ 


r ) 


V 


(<j • t ) 


Y l jm {OA)Xt(i). 


Using the above equation, Eq.(pU|) depends only on the radial coordinates and can be 
expressed as the following integro-differential equation: 
dGu(r) k 


dr . 
dFyjr) 

(IT 

dGyjr ) 

dr N 
dF v (r) 


+ -G v (r) - (E + A - V v (r) + Vs(r))F v (r) +r [ r'dr'A(r, r')F v {r') = 0 
r J 

— —Fu(r) + (E + A — Vv(r ) — Us(r))G t dr) + r / r'dr'Air, r')Gy(r') = 0 
r J 

+ -G v (r) + (E — A + Vy(r) - Vs(r))F v (r) + r [ r'dr'A(r, r')Fu{r') = 0 
r J 


(25) 


-F v (r) - (E - A + VV(r) + U 5 (r))G y (r) + r [ r'dr'A(r,r')Gu(r') = 0, 
dr r J 

where the nucleon mass is included in the scalar potential Vg(r). For the 5-force of Eq.(P^), 

Eq. (^5|) is reduced to normal coupled differential equations and can be solved with shooting 

method by Runge-Kutta algorithms. For the case of Gogny force, the coupled integro- 

differential equations are discretized in the space and solved by the finite element methods. 

The numerical details can be found in Ref.i. Now we have to solve Eqs. (|25|) and ( |TT| ) self- 

consistently for the RCHB case. As the calculation with Gogny force is very time-consuming, 


we solve them only for one case in order to fix the interaction strength for 5-force in Eq.( 


3. The pseudo-spin symmetry 


The Dirac equation in RMF or in the canonical basis of RCHB describes a Dirac spinor with 
mass M moving in a scalar potential Vs(r) and a vector potential Vy(r). With e = M + E, 
the potential V = Vy{r) + Vs(r), which is around —50 MeV, and the effective mass M* = 
M + Vs(r), the relation between the upper and lower components of the wave function can 
be written as: 


1 


9 = 


f = 


E-V 

1 


(d-p) f 


E + 2 M* - V 


(d-p) 9 


(26) 
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Then the coupled equations are reduced to uncoupled ones for the upper and lower com¬ 
ponents, respectively. Effectively we get the corresponding Schrodinger equation for both 
components: 


{a-p)- 


E + 2M--V {3 ' P)9 = {E - V)9 
(3 ■ p)-=2-tf(? ■p)f = (E + 2M* -1 ')f 


(27) 


E — V 

In the spherical case, V depends only on the radius. We choose the phase convention of the 
vector spherical harmonics as: 


(a ■ r)Y l = -Y 

\ ) jm 7 


’ jm 


(28) 


where 


V = 2j-l = 


l + 1, j — l + 1/2 


(29) 


l ~ 1, j - l ~ 1/2 

Here l' is nothing but the pseudo-orbital angular momentum l. After some tedious proce¬ 
dures, one gets the radial equation for the lower and upper components respectively: 


L dr 2 E — V dr dr 

+ [ /v .2 


k (1 — k) 1 hidV lT ^ij 


} F i( r ) 


E — V r dr 
= -(E + 2 M* -V)(E- V)F‘ j (r), 


(30) 


d 2 


d(2M* - V) d 


l dr 2 E + 2 M* — V 

{ n{l + k) 


+ 


dr dr 
k d{2M* - V) 


E + 2 M* — V r dr 

= - (E + 2M* - V)(E - V)G\\r ), 

where 


]G?(r) 


(31) 


k(k — 1 ) = l'{l' + 1 ), k(k + 1 ) = l{l + 1 ). 


(32) 


It is clear that one can use either Eq.(|50|) or equivalently Eq.([H]) to get the eigenvalues E 
and the corresponding eigenfunctions. Normally Eq.(pTl) is used in the literature and the 
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spin-orbital splitting is discussed in connection with the corresponding spin-orbital potential 
1 k d{2M* - V ) 


E + 2 M* — 1 / r dr 

potential ( PSOP ) term, 


If the Eq. 
k dV 


is used instead and the pseudo spin-orbital 


——, is neglected, then the eigenvalues E for the same l 1 


E — V r dr 

will degenerate. This is the phenomenon of pseudo-spin symmetry observed mill. It means 


that Eq.(26) is the transformation between the normal spin formalism and the pseudo-spin 
formalism. 

In Eq. ©, the term which splits the pseudo-spin partners is simply the PSOP. The 
hidden symmetry for the pseudo-spin approximation is revealed as dV/dr = 0, which is 
more general and includes V = 0 discussed in@ as a special case. For exotic nuclei with 
highly diffuse potentials, dV/dr ~ 0 may be a good approximation and then the pseudo-spin 
symmetry will be good. But generally, dV/dr = 0 is not always satisfied in the nuclei and 


the pseudo-spin symmetry is an approximation. However, if 


1 k dV /c(l — k) 


E — V r dr 


•C 


|, the 


pseudo-spin approximation will be good. Thus, the comparison of the relative magnitude of 

W ]_ _ J 

the centrifugal barrier ( CB ),---, and the PSOP can provide us with some information 

on the pseudo-spin symmetry. 

In a recent letteiH, the mechanism behind the pseudo-spin symmetry was studied and 
the pseudo-spin symmetry was shown to be connected with the competition between the 
centrifugal barrier (CB) and the pseudo-spin orbital potential ( PSOP ), which is mainly 
decided by the derivative of the difference between the scalar and vector potentials. With 
the scalar and vector potentials derived from a self-consistent RCHB calculation, the pseudo¬ 
spin symmetry and its energy dependence have been discussed. Here in this paper we will 
extend the previous investigation to the case of exotic nuclei. The pseudo-spin symmetry for 
exotic nuclei is investigated for Zr and Sn isotopes from the proton drip line to the neutron 
drip line. The isospin and energy dependence of the pseudo-spin approximation will be 
investigated in detail in the following section. 
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4. The energy splitting of the pseudo-spin partners 


We use here the non-linear Lagrangian parameter set NLSHo which could provide a good 
description of all nuclei from oxygen to lead. As we study not only the closed shell nuclei, but 
also the open shell nuclei, the inclusion of the pairing is necessary. The pairing interaction 
strength is the same as in Ref.il. The interaction strength in the pairing force of zero range 
Eq.(^) is properly renormalized by the calculation of RCHB with Gogny force. Since we use 
a pairing force of zero range, we have to limit the number of continuum levels by a cut-off 
energy. For each spin-parity channel, 20 radial wave functions are taken into account, which 
corresponds roughly to a cut-off energy of 120 MeV for a fixed box radius R — 20 fm. For 
the fixed cut-off energy and for the box radius A, the strength Vq of the pairing force in 
Eq.(||2j) is determined by adjusting the corresponding pairing energy — |TrAft to that of a 
RCHB-calculation using the finite range part of the Gogny force DISH. We use the nuclear 
matter density 0.152 fm -3 for p 0 . 

The quality of pseudo-spin symmetry can be understood more clearly by considering the 
microscopic structure of the wave functions and the single particle energies in the canonical 
basis. As shown in Ref.!, the particle levels for the bound states in canonical basis are 
the same as those by solving the Dirac equation with the scalar and vector potentials from 
RCHB. Therefore Eqs. ([50]) and ([JI]) are valid in canonical basis after the pairing interaction 
has been taken into account and are very suitable for the discussion of the pseudo-spin 
symmetry. 

The neutron single particle levels in 150 Sn and 120 Zr are given in Fig. la and lb, re¬ 
spectively. The four sets of pseudo-spin partners, i.e., ld 3 / 2 and 2 si/ 2 , I/ 5/2 and 2 p 3 / 2 , 1 ^ 7/2 
and 2d 5 / 2 , 2 d 3 / 2 and 3 si/ 2 , are marked by boxes. As seen in the figure, the energy split¬ 
ting between pseudo-spin partners decreases with the decreasing binding energy. The single 
particle energy of 3 si/ 2 in 120 Zr is —6.00 MeV, and its partner 2d 3 / 2 is —5.86 MeV, the 
splitting is 0.14 MeV. While 2s\/ 2 is —31.62 MeV, 1g? 3 / 2 —33.23 Mev, the splitting is 1.61 
MeV, which is bigger than the former one by a factor of 10. The same situation is found 
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for the energy splitting between pseudo-spin partners in 150 Sn: The single particle energies 

of 3 s!/ 2 and 2d 3 / 2 are —9.645 and —10.11 MeV, respectively. The single particle energies 

for other pseudo-spin partners in 150 Sn are, —11.74 and —13.87 MeV for 2d 5 / 2 and lg-j /2 

partners, —22.46 and —25.50 MeV for 2 p 3 / 2 and I/ 5/2 partners, —33.63 and —36.58 for 2s 3 / 2 

and 1(4/2 partners, respectively. Although we show only the neutron single particle levels 

in 150 Sn and 120 Zr as examples here, the same are found in other Sn and Zr isotopes. It is 

usually seen that the pseudo-spin symmetry approximation becomes better near the Fermi 

surface, which is in agreement with the experimental observation. 

In Fig.l there are also two pairs of pseudo-spin partners ( 3p 3 / 2 and 2 f 5 / 2 partners and 

2 / 7/2 and lhg /2 partners ) near the threshold, apart from the fours pairs of pseudo-spin 

partners below the Fermi level. The energies for these two pairs of pseudo-spin partners 

are —1.581 and —1.031 Mev for the 3 p 3 / 2 and 2 f 5 / 2 partners, and —2.549 and —2.620 MeV 

for the 2 / 7/2 and lhg /2 partners, respectively. Considering their pseudo-spin orbital angular 

momentum 1 — 2 and 4, their splittings A E = 1/ ' 2 ~ - d-f +i /2 are on m —0.1100 

2 / + 1 

MeV and 0.7889 x 10~ 2 MeV, respectively. This is due to the energy dependence and the 

diffuseness of the potential in exotic nuclei, which we will discuss in the following. As it is 

seen in Fig.l, the normal splitting is such that the orbital j = 1 + 1/2 is below the orbital 

j — l — 1/2, except for 3p 3 / 2 and 2 / 5/2 partners. The same also happens for 2 d 3 / 2 and 3 si/2 

partners in Zr isotopes. The pseudo-spin splitting depends on the derivative of the difference 

between the scalar and vector potentials dV/dr, which is small for the exotic nuclei with 

highly diffuse potential. The integration of —^|F| 2 over r gives the splitting of the pseudo- 

dr 

spin partners, whose sign will decide the normal splitting or the reverse. The subtle details 
of the potential are crucial for the pseudo-spin splitting. 

To see the behavior of the pseudo-spin partners around the Fermi level and the isospin 
dependence of the pseudo-spin splitting, we show the single particle levels near the Fermi 
surface in the canonical basis for the Sn and Zr isotopes with an even neutron number as 
a function of the mass number in Fig. 2. The Fermi level is shown by the dashed line. 
The pseudo-spin splitting for the pseudo-spin partners, 2 d 3 / 2 and 3 si/ 2, remains small in 
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Zr and Sn isotopes from the proton drip line to the neutron drip line. The pseudo-spin 
symmetry remains even valid for exotic nuclei. The pseudo-spin symmetry near the neutron 
drip line becomes better than that near the /3-stability line. In Fig. 2a, there is a kink for 
the single particle levels in the continuum, as the contribution from the continuum becomes 
important and the potential becomes diffuse around 130 Sn. But the splitting for 3 p 3 /2 and 
2 / 5/2 partners, and 2 f 7 / 2 and 1 /i 9 / 2 partners in Sn isotopes is small and the pseudo-spin 
symmetry approximation is very good, independent of whether they are in the continuum 
or near the threshold. Therefore we can see that the pseudo-spin symmetry is very well 
reserved for the orbital near the threshod energy and in the continuum region. 


spin orbital splitting more clearly, we plot A E = — ^ 1 ^~ — ^ * +lE versus E = 

for the bound pseudo-spin partners in Sn and Zr isotopes in 


In order to see the energy dependence and the isospin dependence of the pseudo 

n orbital splitting more 

lE Ij=I+ 1/2 + (J + 1/2 

21 + 1 

Fig. 3. In both isotopes, a monotonous decreasing behavior with the decreasing binding 
energy is clearly seen. The pseudo-spin splitting for 3s !/ 2 and 2 d 3 / 2 is more than 10 times 
smaller than that of the 2si/ 2 and ld 3 / 2 . As far as the isospin dependence of the pseudo-spin 
orbital splitting is concerned, the splitting in Sn isotopes gives a monotonous decreasing 
behavior with the increasing isospin. Particularly for 2s 3 /2 and ld 3 / 2 partners, the pseudo¬ 
spin splitting in 170 Sn is only half of that in 96 Sn. Just as we expected, the pseudo-spin 
symmetry in neutron-rich nuclei is better. In Zr isotopes, although the situation is more 
complicated ( e.g., the effect of the deformation which is neglected here ) , the pattern is 
more or less the same, i.e., a monotonous decreasing behavior with the decreasing binding 
energy and a monotonous decreasing behavior with the isospin. From these studies, we see 
that the pseudo-spin symmetry remains a good approximation for both stable and exotic 
nuclei. A better pseudo-spin symmetry can be expected for the orbital near the threshold, 
particularly for nuclei near the particle drip-line. 
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5. The pseudo-spin orbital potential 


To understand why the energy splitting of the pseudo-spin partner changes with different 

binding energies and why the pseudo-spin approximation is good in RMF, the PSOP and 

CB should be examined carefully. Unfortunately, it is very hard to compare them clearly, as 

the PSOP has a singularity at E ~ V. As we are only interested in the relative magnitude 

of the CB and the PSOP, we introduce the effective CB, (E — V )---, and the effective 

k dV 

PSOP, ——, for comparison. They correspond to the CB and the PSOP multiplied by a 
r dr 

common factor E — V respectively. 

The effective PSOP does not depend on the binding energy of the single particle level, 
but depends on the angular momentum and parity. On the other hand the effective CB 
depends on the energy. Comparing these two effective potentials one could see the energy 
dependence of the pseudo-spin symmetry. They are given in Fig. 4 for s \/2 ( lower ) and 
d' 3/2 ( upper ) of 120 Zr in arbitrary scale. 

The pseudo-spin approximation is much better for the less bound pseudo-spin partners, 
because the effective CB is smaller for the more deeply bound states. This is in agreement 
with the results shown in Fig.3. The effective PSOP and the effective CB are also given as 
inserts in Fig.4 in order to show their behavior near the nuclear surface. 

In order to examine this carefully, we compare the effective CB ( dashed lines or dot- 
dashed lines ) and the effective PSOP ( solid lines ) multiplied by the squares of the lower 
component wave function F(r), which are given in Fig. 5, for 2 si/ 2 ( upper left ), 3 si/ 2 ( 
lower left), ld $/2 ( upper right ), and 2 ^ 3/2 ( lower right ) of 120 Zr in arbitrary scale. The 
pseudo-spin approximation is much better for the less bound pseudo-spin partners, because 
the effective CB is smaller for the more deeply bound states. This is in agreement with the 
results shown above. The integrated values of the potentials in Fig.2 with r are proportional 
to their contribution to the energy after some proper renormalization. It is clear that the 
contribution of the effective CB ( dashed lines or dot-dashed lines ) is much bigger than 
that of the effective PSOP ( solid lines ). Generally the effective PSOP is two orders of 
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magnitude smaller than the effective CB. 

In Fig. 1 and 2 we notice that the orbital j = l + 1/2 is generally below the orbital 

j = l — 1/2, except for 3 p 3 / 2 and 2/ 5 / 2 partners in Sn isotopes. The same situation also 

happens for 2 d 3 / 2 and 3 si /2 partners in Zr isotopes. As the pseudo-spin splitting depends 

on PSOP, which depends on the subtle radial dependence of the potentials, sometimes the 

PSOP may have positive or negative regions as a function of r which cancell each other. The 

integration of ——I FI 2 over r gives the splitting of the pseudo-spin partners, whose sign will 
dr 

decide the normal splitting or the reverse. That is the reason why the orbital j = l + 1/2 is 
above the orbital j = l — 1/2 for 3p 3 /2 and 2 / 5/2 partners in Sn isotopes and for 2d 3 / 2 and 
3 si /2 partners in Zr isotopes. 

6. The wave-function of pseudo-spin partners 

In the above discussion, we have seen that the PSOP is much smaller than the CB. Therefore 
if we neglect the PSOP in Eq.(^D|), the lower component of the Dirac wave functions for the 
pseudo-spin partners will be the same, i.e., in the case of the exact pseudo-spin symmetry, 
the lower component of the pseudo-spin partners should be identical (except for the phase). 
The upper component of the Dirac wave functions can be obtained from the transformation 
in Eq. (j26[) , which depends on the quantum number k. Therefore the study of the Dirac wave 
functions for the pseudo-spin partners will provide a check for the pseudo-spin approximation 
in nuclei. As examples, the normalized single nucleon wave functions for the upper ( G ) and 
lower ( F ) components of the Dirac wave functions for the pseudo-spin partners ld 3 / 2 and 
2 si/ 2 , I/ 5/2 and 2p 3 / 2 , 1 ^ 7/2 an d 2d 5 / 2 , 2d 3 / 2 and 3 si / 2 in 120 Zr are given in Fig. 6 . Of course, 
the lower components are much smaller in magnitude compared with the upper component 
in Eq.(^6|). The phase of the Dirac wave functions for one of the pseudo-spin partners has 
been reversed in order to have a careful comparison. It is seen that the lower components of 
the Dirac wave functions for the pseudo-spin partners are very similar and are almost equal 
in magnitude, as observed also for 208 Pb in Refill. The similarity in the lower components 
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F of the wave function for the pseudo-spin partners near the Fermi surface is better than 
for the deeply bound ones. The lower components for the pseudo-spin partners with small 
pseudo-spin orbital angular momentum are better than for the ones with large pseudo-spin 
orbital angular momentum. As seen in Fig. 6 , the similarity for pseudo-spin partners 2 d 3 /2 
and 3 si / 2 is better than for pseudo-spin partners ld 3 / 2 and 2 si/ 2 . The similarities for pseudo¬ 
spin partners, 2 d 3 / 2 and 3 si/ 2 , 1 d 3 / 2 and 2 si/ 2 , are better than for the pseudo-spin partners 


I/ 5/2 and 2p 3/2 , lg 7/2 and 2 d 5/2 . 

Although the lower components for the pseudo-spin partners are very close to each other, 
the difference for the upper components is very big. The upper component of the Dirac wave 
functions can be obtained from the transformation in Eq. (|26| ) , which for the sperical case 

can be reduced to the follows: 

1 r dF l J(r) k 


dh r) = 


+ 

r 


E-V L dr ' r~ l v ' J ' ^ 

As seen in Fig.6, in the case of exact pseudo-spin symmetry, where both E and F l d (r) are 

identical for the pseudo-spin partners, the upper conpnonents G 1 / (r) will be different due to 

the term —F l Mr). For the pseudo-spin partners with small l, the contribution of the term 
r 

— F l Mr ) becomes less important for larger r and a similarity between the upper components 
r 

can happen in the nuclear surface. As for examples, for r > 6 fm, the upper components 
for the pseudo-spin partners ld 3 / 2 and 2 si/ 2 , I/ 5/2 and 2 p 3 / 2 , 2 d 3 / 2 and 3 si/ 2 , in Fig. 6 are 
very similar. 


7. Summary 

In conclusion, the pseudo-spin symmetry is examined in normal and exotic nuclei in the 
framework of RCHB theory. Based on RCHB theory the pseudo-spin approximation in exotic 
nuclei is investigated in Zr and Sn isotopes from the proton drip line to the neutron drip line. 
The quality of the pseudo-spin approximation is shown to be connected with the competition 
between the centrifugal barrier (CB) and the pseudo-spin orbital potential ( PSOP ), which 
is mainly decided by the derivative of the difference between the scalar and vector potentials 
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dV/dr. If the derivative of the difference between the scalar and vector potentials dV/dr 
vanishes, the pseudo-spin symmetry is exact. The condition dV/dr ~ 0 may be a good 


approximation for the exotic nuclei with highly diffuse potential. Further the new condition 
1 k dV k(1 — k) 


E — V r dr 




is found under which the symmetry is preserved approximately. 


We have examined this condition to see how good the pseudo-spin symmetry is in RCHB. For 

a given angular momentum and parity channel, the effective CB, (E — V) ---, becomes 

stronger for the less bound level, so the pseudo-spin symmetry for the weakly bound state 
is better than that for the deeply bound state, which is in agreement with the experimental 
observation!®. The pseudo-spin symmetry is found to be a good approximation even for the 
exotic nuclei with highly diffuse potential. The above conclusion has been well supported by 
RCHB calculations for Zr and Sn isotopes from the proton drip line to the neutron drip line. 
From the simple Dirac equation, it has been shown that there are two equivalent ways to 
solve the coupled Dirac equation for the upper and lower components: i.e., the normal spin 
formalism and pseudo-spin formalism. Both formalisms are equivalent as far as the energies 
and wave functions are concerned. Their relation is given by Eq. (|2E|) , which indicates that 
the unitary transformation from the conventional formalism to the pseudo-spin formalism 
has the ”p-helicity”000. Summarizing our investigation, we conclude: 


1 The quality of the pseudo-spin approximation is connected with the competition be¬ 
tween the CB, and the PSOP which is mainly proportional to the derivative of the 
difference between the scalar and vector potentials dV/dr ; 

2 The pseudo-spin symmetry is a good approximation for normal nuclei and become 
much better for exotic nuclei with highly diffuse potentials; 

3 The pseudo-spin symmetry has strong energy dependence. The energy splitting be¬ 
tween the pseudo-spin partners is smaller for orbitals near the Fermi surface. 

4 The energy difference between the orbital j — l + 1/2 and the orbital j = l — 1/2 is 
always negative, except for 3p 3 / 2 and 2 f 5 / 2 partners. The same situation also happens 
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for 2ds/2 and 3 si/ 2 partners in Zr isotopes. The integration of -j—\F\ 2 over r gives the 
splitting of the pseudo-spin partners, whose sign will decide the normal splitting or 
the reverse. 


5 The lower components of the Dirac wave functions for the pseudo-spin partners are 
very similar and almost equal in magnitude. The similarity in the lower components 
of the wave function for the pseudo-spin partners near the Fermi surface is closer than 
for the deeply bound ones. 
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Figure Captions 


Fig. 1 The single particle levels in the canonical basis for the neutron in 120 Zr and 150 Sn. 
The Fermi surface is shown by a dashed line. The bound pseudo-spin partners are 
marked by boxes 


Fig. 2 The single particle energies of the neutron in the canonical basis as a function of the 
mass number for Zr and Sn isotopes. The dashed line indicates the chemical potential. 


Fig. 3 The pseudo-spin orbit splitting A E = 


ergy E = 


lEij=i+ 1/2 + (^+ l ) E ij=i -\/2 


E lj=l- 1/2 E lj=l+ 1/2 

21 + 1 


versus the binding en- 


21 + 1 


for Zr and Sn isotopes. From left to right, 


the pseudo-spin partners correspond to (ld 3 / 2 , 2 si/ 2 ), (I/5/2, 2p 3 / 2 ), (1^7/2, 2d 5 / 2 ) and 
(2d 3 / 2 , 3 si/ 2 ), respectively. 


Fig. 4 The comparison of the effective centrifugal barrier ( CB ) (E — V) 


k,(k — 1 ) 


dashed 


lines and dot-dashed lines ) and the effective pseudo-spin orbital potential ( PSOP ) 

k dV 

—— ( solid line ) in arbitrary scale for d 3 / 2 ( upper ) and Si / 2 ( lower ) in 120 Zr. The 

r dr 

dashed lines are for ld 3 / 2 and 2 si/ 2 , and the dot-dashed lines are for 2d 3 / 2 and 3 si/ 2 . 
The inserted boxes show the same quantities, but the ordinate is magnified and the 
abscissa is reduced to show the behaviors of the effective CB and the effective PSOP 
near the nuclear surface. 


Fig. 5 The comparison of the effective centrifugal barrier ( CB ) {E — V) 


— 1 ) 


dashed 


lines and dot-dashed lines ) and the effective pseudo-spin orbital potential ( PSOP 

n dV 

)-( solid line ) multiplied by the square of the wave function F of the lower 

r dr 

components in arbitrary scale for d 3 / 2 ( upper ) and Si / 2 ( lower ) in 120 Zr. The dashed 
lines are for ld 3 / 2 and 2 si/ 2 , and the dot-dashed lines are for 2d 3 / 2 and 3 si/ 2 . 


Fig. 6 The upper component G and lower component F of the Dirac wave functions for 
the pseudo-spin partners in 120 Zr. The phase of the Dirac wave functions for one of 
the pseudo-spin partners has been reversed in order to have a careful comparison. 
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